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ABSTRACT 

We use three-dimensional hydrodynamic simulations to investigate the nonlinear gravitational re- 
sponses of gas to, and the resulting drag forces on, very massive perturbers moving on circular orbits. 
This work extends our previous studies that explored the cases of low-mass perturbers on circular 
orbits and massive perturbers on straight-line trajectories. The background medium is assumed to be 
non-rotating, adiabatic with index 5/3, and uniform with density po and sound speed ag. We model 
the gravitating perturber using a Plummer sphere with mass M p and softening radius r s in a uniform 
circular motion at speed V p and orbital radius R p , and run various models with differing 1Z = r s /R p , 
M = V p /ao, and B = GM p / (a^Rp). A quasi-steady density wake of a supersonic model consists of a 
hydrostatic envelope surrounding the perturber, an upstream bow shock, and a trailing low-density 
region. The continuous change in the direction of the perturber motion makes the detached shock 
distance reduced compared to the linear-trajectory cases, while the orbit-averaged gravity of the per- 
turber gathers the gas toward the center of the orbit, modifying the background preshock density to 
p 1 k, (1 + 0.46i3 1 ' 1 )po depending weakly on A4. For sufficiently massive perturbers, the presence of 
a hydrostatic envelope makes the drag force smaller than the prediction of the linear perturbation 
theory, resulting in F = Anpi(GM p ) 2 /V 2 x (0.7rj~ r ) for rj B = B/(M 2 - 1) > 0.1; the drag force for 
low- mass perturbers with r\ e < 0.1 agrees well with the linear prediction. The nonlinear drag force 
becomes independent of 1Z as long as 1Z < ?7e/2, which places an upper limit on the perturber size for 
accurate evaluation of the drag force in numerical simulations. 
Subject headings: hydrodynamics — ISM: general — shock waves 
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1. INTRODUCTION 

Dynamical friction (DF) arising from the gravitational 
interaction of a massive body with a background medium 
is considered to be one of the most powerful mechanisms 
respo nsible for the orbital decay of astronomical objects 
(e.g., iBinnev fc Tremaind 120081 and references therein). 
DF occurs not only in a collisionless background in which 
long-range, two-body interactions transfer momentum 
from an object in motion to the field particles, but also 
in a continuous gaseous medium where a perturber loses 
its momentum due to the gravitational drag force ex- 
erted by its own induced wake. The DF in a gaseous 
background is likely to play an important role in the 
orbital decay of companions in common-envelope bina- 
ries, supermassive black holes (SMBHs) at galaxy cen- 
ters, massive galaxies in galaxy clusters, etc. For in- 
stance, in a binary system, DF causes a low-mass com- 
panion engulfed by an extended primary to spiral in to- 
ward the core of the primary, which in turn spins up 
the primary envelope, greatly affecting the subsequent 
evolution of the system (e.g.. iTaam fe S andquiste I200C 



Edgar! 12004 iNordhaus fc Blackmanl 120061: iMaxted et al.l 
20091 and references therein). Also, numerical simulations 
show that the DF due to a gaseous medium expedites 
the g rowth of SMBHs by mergers i n colliding galaxies 
(e.g.. iMilosavlievic fc Merrit l] |2003j: lEscala et all [20M - 
2005 [lDotti et al.ll2006l : iMaver et al.l l2007: Colpi & Dotti 
20091) . potentially explaining the ubiquity of SMBHs at 
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galactic nuclei (e.g.. iBegehnan et al"1ll980l : iMenou et al.1 
l200lt iFerrarese fc Fordll2005ri . 

Using a tim e-dependent linear perturbation theory, 
lOstrikerl t|1999f ) showed that the drag force on a perturber 
with mass M p moving at speed V p on a rectilinear tra- 
jectory through a uniform gaseous medium with density 
Po and sound speed ap is given by 



Fr, 



4ttpq(GM p ) 2 

V 2 
v 



±ln 



l+M 
l-M 



-M, 



where Ai = V p /a Q is the Mach number, and 



M<1, 

M>1, 

(1) 
is the 



characteristic size of the perturber. The key features 
of equation (fTJ) are that (i) the gaseo us DF force be- 
comes identical to the classical formula of lChandrasekharl 
(|1943|) for the collisionless drag when M. 3> 1, (h) the 
gaseous DF force is more efficient than the collisionless 
drag for M. ~ 1, and (iii) the gaseous DF force is non- 
vanis hing even for subso nic perturbers with M. < 1 (see 
also Uust fc Kegellll990D . The last point improves the 
previous notion that the gaseous drag is absent in sub- 
sonic cases due to the front-back symmetry in the steady- 
state wakes (e..g.lDokuchaevlll964l iRuderman fc Spiegell 
119711: iRephaeli fc Salpeterlll980l) . Several recent studies 
showed that equation (fTJ) can be applicable, with some 
modifications, to more general cases such as in a radi- 
ally s tratified medium ([Sanchez-Salcedo fc Bra ndenburg 
2001tl. for circu lar-orbit perturbers (|Kim fc Kiml 120071 : 



Kim et all 
(Bar ausseT 



2008), for perturbers with relativi stic speed 
2007) or in accelerating motion (Namouni 
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|2010|) . etc. In particular. iKim fc Kiml (|2007l hereafter 
KK07) used a semi-analytic method to show that equa- 
tion (jj is a good approximation to the DF force on 
circular-orbit perturbers provided V p t = 2R p , where R p 
is the orbital radius. 

While all the theoretical studies mentioned above con- 
sider low-mass pertu rbers and find various astrophysi 
cal applications ( e .g.. iNaravan! 2000HE1-Zant et all '. 
Kim et all 120051: IKiml 120071: IConrov fc Ostrikerl E 
Villaver fc Livioll2009( ). there are some situations such as 



2008 



in orbital decay of SMBHs or companions in common- 
envelope binaries, where perturbers have so large masses 
that the induced density wakes are in t he nonlinea r 
regim e. Using hydrodynamic simulations, IKim fc Kir 



Kim 
rikerl 



2009, hereafter KK09) extended the work of lOstriker 
1999) to study nonlinear DF force for a very massive 
perturber on a straight-line trajectory. By modeling a 
perturber using a Plummer sphere with softening radius 
r s , KK09 found that a nonlinear supersonic wake is char- 
acterized by a detached bow shock and a hydrostatic en- 
velope near the perturber. The resulting drag force de- 
pends solely on the dimensionless parameter r] A defined 
as 



Va = 

and is given by 
F 

F~ 



A 



M 2 



1 



with A 



GM n 



-0.45 



for j] A > 2, 



(2) 



(3) 



while F/Fu n 1 for r\ A < 0.7. The reduction of the 
nonlinear DF force compared to the linear estimate is 
due to the presence of a hydrostatic envelope that makes 
the density distribution spherically symmetric near the 
perturber. This clearly demonstrates that the nonlinear 
effect can be significant if a perturber is very massive. 

Since astronomical objects usually follow curvilinear 
rather than straight-line trajectories, it is interesting to 
see whether equation ([3]) remains valid for perturbers on 
circular orbits. For this purpose, we in this paper take 
one step further from KK09 to study the nonlinear gravi- 
tational responses of gas to, and the associated drag force 
on, a very massive perturber moving on a circular orbit. 
This work also extends the linear perturbation analyses 
of KK07 by considering the nonlinear effect on the den- 
sity wakes. The intensity of gravitational influence of a 
perturber to the background gas can be measured by the 
Bondi radius re = GMp/a 2 ]. For perturbers on straight- 
line trajectories, the softening radius (or the perturber 
size) is the lone length scale with which the Bondi ra- 
dius can be compared. The drag force, correspondingly, 
depends on r s and M p only through the dimensionless 
parameter A: increasing r s is equivalent to decreasing 
M p . Equation ([3]) then predicts that the drag force on a 
highly nonlinear object with A 3> 1 would be negligible, 
which results simply from the fact that r s is inseparable 
from M p . On the other hand, circular orbits with orbital 
radius R p naturally introduce an additional length scale, 
so that the dependence of the drag force on r s can be 
explored independently of the perturber mass. 

While a perturber in real astronomical situations is 
likely to move through a background medium that is 
rotating and/or stratified, we in this work idealize the 
gaseous medium as being initially static and uniform. 



This certainly introduces a few important caveats that 
should be noted from the outset. If the gaseous medium 
is supported primarily by rotation, as in protoplane- 
tary disks, the perturber launches spiral waves at Lind- 
blad resonances that propagate away from it (e.g.JWardl 
19971 : iTanaka et al.ll2002t ICIiambersl[2009t iLubow fc Idal 
20101 and references therein), which is quite different 



from quasi-steady density wakes with long trailing tails 
produced in a non-rotating medium (e.g., KK07). If 
the medium is instead supported by thermal pressure, 
it should have a stratified density in the radial direc- 
tion, as in intracluster media or common-envelope bi- 
naries. In this case, the gradient in the background 
density profile can be ignored if the Bondi radius is 
much smaller than both the pressure scale length and 
the orbital radius, that is, if M < (M(R p )/M p ) 1 ^ 
and B < 1, where M(R p ) is the dynamical mass en- 
closed within R p and B is the dimensionless perturber 
mass defined in Sj2] below. Although the first condition 
is readily met, for example, in common-envelope bina- 
ries with planets or brown dwarfs as lo w-mass compan- 
ions (e.g.. lNordhaus fc Blackmanll2006t ). the second con- 
dition is not well satisfied especially for very massive per- 
turbers with B > 1 considered in this paper. In the latter 
case, the wakes and the associated DF forces are likely 
to be affected by t he density gradient of the background 
medium (e.g., see Dust fc Penarrubiall2005l for collision- 
less cases). Neglecting the density stratification in the 
background medium also suppresses gas buoyancy, pre- 
cluding the potential effects of convective motions and 
gravity modes in heating the medium by resonant ex- 
citations (e.g.. iBalbus fc Sokerl H990: Luf kin et all 119951: 
lKimll2007l) and other processes. 

In addition to the above assumptions, we ignore the 
orbital motion of the background gas with respect to 
the center of mass of the whole system, which can be 
a valid approximation only when B/M 2 < 1. If this 
condition is not satisfied, the neglect of the centrifugal 
and Coriolis forces arising from the orbital motion of the 
background may affect the density wakes and the drag 
forces (e.g.. lAdams et al.lfl989l: lOstriker et al.lll992| ). We 
also treat the gas using an adiabatic equation of state, 
which implicitly assumes that the orbital energy of the 
perturber is converted to heat as it spirals inward, and 
thus is valid only if radiative loss is negligible. By ignor- 
ing self-gravity, we do not consider any back reaction of 
the gas on the perturber. 

Given these limitations and constraints, we by no 
means attempt to apply the results of this work to real 
astronomical systems. Nevertheless, the idealized mod- 
els considered in this paper help to isolate the effect 
of the perturber size or orbital radius on the nonlin- 
ear DF force. The results of this work will be partic- 
ularly useful to justify large perturber sizes employed in 
recent hydrodynami c simulations, such as for SMBHs at 
galac t ic nuclei (e.g..lEscala et al. 2004, 2005; Dotti et all 
120061 [20071 iMaver et all 120071 : iCuadra et all l2009irand 
companions in common-envelope binaries (e.g.. [ Ruffertl 
119931 : ISandquist et all Il998l : iRicker fc Taamll2008l ). etc. 
These simulations usually treat the perturber using a 
softened point mass, with its size inevitably limited by 
numerical resolution. For instance, iV-body/SPH simu- 
lations for the orbital decay of SMBHs take quite large 
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values (up to a fe w pc) for r s (e.g., lEscala et al.l 12004 
iMayer et al.ll2007f ) , although the realistic values are prob- 
ably of order of the Schwartzschild radius (~ 1 AU 
for M p — 2 x 10 7 Mq). Grid-based simulations of the 
common-envelope phase of binaries also represent a com- 
panion using a point mass with size set by the grid spac- 
ing, which is larger than the actual co mpanion size by 
one or two orders o f mag nitude (e.g., ISandauist et al.l 
119981: iRicker fe Taam1l2008t) . If the drag force on circular- 
orbit perturbers depends on r s similarly to the linear- 
trajectory cases, simulations with such large r s would 
significantly underestimate the real decay time since the 
induced density wake would erroneously be in the linear 
regime. If the drag force instead turns out insensitive to 
r s in circular-orbit cases, a large value of r s would be 
reasonable as long as it does not change the drag force 
much. 

The paper is organized as follows: In fj2] we describe 
numerical methods and models we adopt. In fj3j we run 
models for low-mass perturbers and compare the result- 
ing drag forces with the analytical predictions to find 
a proper relationship between r s and r m - ln . In ij4] we 
present evolution and quasi-steady distributions of fully 
nonlinear density wakes and the associated drag forces 
on massive, circular-orbit perturbers. Finally in ij5l we 
summarize our findings and discuss their implications on 
the proper choice of the perturber size. 

2. NUMERICAL METHOD 

We consider a massive perturber moving on a circular 
orbit through a gaseous medium, and study the gravita- 
tional response of the gas to the perturber and the result- 
ing drag force. Similar work for low-mass perturbers on 
circular orbits and massive perturbers on straight-line 
trajectories was presented in KK07 and KK09, respec- 
tively. Since the DF timescale is usually much longer 
than the orbital time, the idealized circular orbit is 
a reasonable approximation to real curvilinear orbits. 
We assume that the background medium is unmagne- 
tized, non-self-gravitating, and initially static and ho- 
mogeneous with density po and adiabatic sound speed 
ao. We adopt an adiabatic equation of state with index 
7 = 5/3. The perturber with mass M p moves with a 
constant angular velocity f2 p along a circle with radius 
R p on the z = plane; the corresponding linear veloc- 
ity is V p = R p ft p - We represent the perturber using a 
Plummer potential 



<f> P (x,t) = — 



GM„ 



>(*)! 



for t > 0, (4) 



where r s is the softening radius and x p (t) — 
Rp(cos£lpt,smQ p t,0) is the perturber location at time 
t in the Cartesian coordinates (x,y,z). 

We take R p , ao, and t cross = R p /ao as the units 
of length, velocity, and time, respectively, in our sim- 
ulations. Then, our models are completely character- 
ized by three dimensionless parameters: 1Z = r s /R p , 
M = V p /a , and 



B _ GMp 



TLA, 



(5) 



which is the ratio of the Bondi radius to the orbital ra- 
dius. Note that 1Z is a new parameter introduced by the 



circular orbit, which is absent for perturbers on straight- 
line trajectories. For future purposes, we define a dimen- 
sionless parameter 

similarly to r/ A in equation @. To explore the paramet- 
ric dependence of the drag force on K, M, and B (or 
equivalently A), we run a total of 65 three-dimensional 
simulations with 1Z varying in the range of 0.025-0.4, M. 
in 0.5-4, and B from 5 x 10~ 4 to 4 (A from 0.01 to 80). 
Our fiducial model has K = 0.05, M = 2, and B = 0.5. 
We use the orbital time £ or b = 2ttR p /V p — 2nt CIOSS /M 
as the time unit of our presentation. 

We integrate the basic equations of ideal hydrody- 
namics using a modifi ed version of the ZEUS code 
(|Stone fe Normanl Il992f), parallelized on a distributed- 
memory platform. Our simulations domain is a rectangu- 
lar box which spans the region — 10i? p < x,y, z < 10R p . 
We construct a logarithmically spaced Cartesian grid 
with 512 3 zones, with spacings of 0.0075i? p , 0.018i? p , and 
0.11i? p at \x\/R p = 0, 1, and 10, respectively. For some 
selected parameters, we have also run low-resolution sim- 
ulations with 256 3 zones and checked that the resulting 
drag forces agree, within a few percents, with those from 
the high-resolution models. We adopt the outflow bound- 
ary conditions at all the boundaries and run the simula- 
tions typically up to 15i croS s- The simulation outcomes 
are insensitive to the box size and the boundary condi- 
tions we adopt since it takes about ~ 20i cross for sound 
waves to travel across the simulation box, while the drag 
forces saturate typically within ~ 0.5i or b ~ (2 — 3)t croS s- 

3. LINEAR CASES 

The time-dependent linear-perturbation analyses for 
the wakes of low-mass perturbers usually consider a 
point- mass object with vanishing r s in equation (j4j), 
which in turn requires to introduce the cut-off radius 
r m j n i n the linear D F force formula (e.g, Uust fe Kegell 
119901: lOstrikerl [l99l KK07). On the other hand, one 
has to adopt a non-zero value for the softening radius 
in numerical simulations, making the DF force depen- 
dent on r s in the linear regime. To explore how the per- 
turber mass affects the drag force in comparison with the 
linear theory, therefore, it is necessary to find a proper 
relationship between r s and r m i n that makes the numeri- 
cal results consistent with the analytic predictions when 
B <C 1. For this purpose, we present in this section the 
results of numerical simulations for low-mass perturbers 
with B = 5x 10~ 4 (A = 0.01) and K = 0.05 as functions 
otM. 

Figure [T] shows the snapshots of density wakes on the 
x-y plane at f/i OT b — 1 for M = 0.8, 1.5, 2, and 3 cases. 
The perturber in each frame is moving along the white 
circle in the counterclockwise direction. For subsonic 
perturbers, the circular orbit makes the density wakes 
curved along their orbits that would otherwise remain 
symmetric with respect to the line of motion. For super- 
sonic models, the perturber is able to overtake its own 
wake, creating a strong trailing tail bounded by Mach 
waves (KK07). As the Mach number increases from 
unity, the opening angle of the head of the curved Mach 
cone deceases, while the tail thickens. The dotted lines 
for the supersonic models plot equation (B2) of KK07 
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Fig. 1. — Steady-state distributions of the perturbed density in logarithmic color scale on the x-y plane at t/t rb = 1 f° r low-mass 
perturbers with 1Z = 0.05 and B = 5 X 10 -4 (A = 0.01). Colorbar labels log(p/po — !)• The white circle in each panel denotes the orbit, 
while the small black circle near (x,y) = (R p ,0) marks the softening radius of the perturber. The spiral dotted lines for the supersonic 
cases draw the outer edges of the tails from the linear theory, eq. (B2) of KK07, which are in good agreement with the numerical results. 



for the shape of the tail 1 , which are in excellent agree- 
ment with the outer boundaries of the tails produced in 
the simulations. Nevertheless, the finite perturber size in 
the latter makes the tail boundaries broader compared 
to the point-mass counterparts with r s — 0. 

For a given density wake p{x 1 t), it is straightforward 
to calculate the gravitational drag force on the perturber 
by evaluating the integral 

_ f GM P (\x - x p \)p(x,t)(x - x p ) ■ 3 
1 ' ~J \x-x p \z ' 1 ' 

where M P {r) = M p r 3 /(r 2 + r 2 ) 3 / 2 is the mass distribu- 
tion of the Plummer sphere and tp = z x x^ is the unit 
vector in the azimuthal direction. Figure [zp shows the 

1 Equation (B2) in KK07 reduces to R/R p = M~ 1 s for R/R p > 
1, with s being the azimuthal angle, suggesting that the wake tail 
at large radii resembles an Archimedean spiral. 



temporal evolution of F(t) for the cases shown in Fig- 
ure [T] The friction force saturates to a constant value 
in less than ~ 2t cross ; the amplitudes of fluctuations in 
F are less than 1% of the mean values. This is unlike 
in the linear-trajectory models where the DF force on a 
supersonic perturber increases logarithmically with time, 
since the whole density wake growing in size with time is 
located behind the perturber. The region of influence 
in the circular-orbit models expands at a sonic speed 
from the orbit center. Consequently, the far-field wake in 
these models is more or less spherically symmetric and 
thus does not contribute to the DF force much. Fig- 
ure^ plots as solid circles the steady-state drag forces 
against M. for low-mass perturbers with B = 5 x 10~ 4 
(A = 0.01) and 1Z = 0.05. The solid line represents the 
fit using equation ([l} with 

r mi n/r s = 1.5 + 0.7(X-1.8) 2 , (8) 
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Fig. 2.— (a) Temporal evolution of the dimensionless DF force for TZ = 0.02, B = 5 X 10~ 4 (.4 = 0.01), and M = 0.8, 1.5, 2.0, and 3.0. 
The drag force reaches a steady-state value within ~ 2t cr oss. (&) The mean drag force (solid circles), averaged over t/t CTOBB = 2-5, against 
the Mach number. The standard deviation is smaller than the size of the circles. The solid line plots the linear force formula (eq. [T]) with 
V p t = 2R P and r m i n /r s = 1.5 + 0.7(M — 1.8) 2 in good agreement with the numerical results. 



which is in good agreement with the numerical results. 
In what follows, we will use equation ([5]) as the cut-off 
radius when we compare the numerical DF forces with 
the analytic estimates. 

4. NONLINEAR CASES 

4.1. Wake Evolution 

We now focus on the cases with very massive per- 
turbers. Figure [3] plots the snapshots of logarithmic 
density on the x— y plane for our standard model with 
Tl = 0.05, M = 2, and B = 0.5 (A = 10). In each 
panel, the white line draws the perturber orbit, while the 
small circle in black marks the perturber location with 
its size corresponding to the softening radius r s . The 
perturber is initially located at (x,y) = (R p ,0) and or- 
bits in the counterclockwise direction. Note that the size 
of the region displayed in Figure [3] differs from panel to 
panel for clear presentation. The evolution of the den- 
sity wake at early phase is qualitatively similar to the 
linear-trajectory cases described in KK09. An introduc- 
tion of the perturber at t = provides the background 
gas with strong perturbations that readily develop into 
a bow shock. The incident gas along the line of motion 
is shocked and gathered near the perturber. This sets 
up a strong pressure gradient in the region between the 
shock and the perturber, which in turn causes the shock 
to slowly move away from the perturber in the upstream 
direction until it reaches an equilibrium position. Some 
gas flowing with a non-zero impact parameter is grad- 
ually pulled by the perturber, adding to the density at 
the rear side (Fig. Due to the strong gravity, the 

material piled up at the backside is pulled back toward 
the perturber, creating a strong counterstream along the 
line of motion as well as a surrounding vortex ring (see 
KK09). 



For the linear-trajectory cases, KK09 has shown that 
the counterstream slows down due in part to the strong 
pressure gradient established near the perturber and in 
part to a buoyant expansion of the vortex ring in the 
lateral direction. The interaction of the counterstream 
with the shocked gas causes the bow shock to move back 
and forth around an equilibrium position. During this 
time, the vortex ring also undergoes an overstable oscil- 
lation whose amplitude grows secularly with time. When 
the vortex ring moves out beyond a half of the Hoylc- 
Lyttleton radius thl = 2GM p /Vp, it becomes less gravi- 
tationally bound and soon swept downstream by the ram 
pressure of the incident gas, leaving a near-hydrostatic 
envelope surrounding the perturber. This is how a den- 
sity wake enters into a quasi-steady state in the linear- 
trajectory cases. For the circular-orbit cases, however, 
the line of perturber motion keeps changing with time. 
With the perturber displaced from the direction of the 
counterstream, the counterstream is virtually unhindered 
and travels almost straight in the direction tangent to 
the orbit, producing the protruding region in the wake 
at x ~ y ~ 0.85R P shown in Figure [3l> . Consequently, 
the vortex ring manifesting its presence by low density 
at x/R p = 0.6-0.9 and y/R p = 1.1-1.2 in Figured that 
was carried by the counterstream decouples from the per- 
turber, promptly leaving a quasi-static density distribu- 
tion around the perturber. Note that the time to form a 
near-hydrostatic envelope is about ~ 0.5i or b ~ 30r s /ao, 
which is shorter than in the linear-trajectory cases by 
about an order of magnitude. 

While the density distribution near the perturber is 
nearly spherical, a part of the wake outside the orbit is 
pushed radially outward by the counterstream, develop- 
ing a trailing tail on the x— y plane. As the counterstream 
weakens, the low-density regions associated with the vor- 
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Fig. 3. — Density snapshots on the x—y plane for a supersonic model with M. = 2, 1i = 0.05, and B = 0.5 (.4 = 10). The size of the 
region displayed differs from panel to panel for clarity. The white circle or part of it in each panel denotes the orbit, while the small black 
circle marks the softening radius of the perturber. In (/), the black dotted line plots equation (B2) of KK07 for the linear tail, while the 
white dashed line that matches the outer edge of the nonlinear tail is drawn by rotating the former by 75° in the counterclockwise direction. 
Colorbar labels log(p/po — !)■ See text for detail. 
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Fig. 4. — Quasi-steady distributions of the density wake and velocity field, seen in the inertial frame, on the x/R p = 1 (bottom right), 
y = (top left), and z = (bottom left) planes for a supersonic model with M = 2, 1Z = 0.05, and B = 0.5 (A = 10) at t/*orb = 2. The 
small black circle marks the softening radius of the perturber, while the white curve in the bottom left panel draws the perturber orbit. 
The length of the line segment shown in the upper left corner of each panel, corresponding to twice the sound speed, measures the size of 
the velocity vectors. Colorbar labels log(p/po — 1). 
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Fig. 5. — Distributions of quasi-steady density wakes with differing J\A on the x—y plane, with the perturber located at (x, y, z) = (R p , 0, 0). 
All the models have 1Z = 0.05 and 8 = 1. The white circle in each panel shows the orbit. Colorbar labels log(p/po — !)• 



tex ring first merge together and are then stretched along 
the orbit as the perturber continues the orbital motion 
(Fig. e )- In the low- mass perturber cases with A <C 1, 
the wake tail is created by the overlapping of the Mach 
cone and the sonic sphere (KK07). In the nonlinear cases, 
however, it is the low-density regions produced by the 
counterstream that separate the tail from the rest of the 
wake. The black dotted curve in Figure [3)f plots the 
shape of the linear tail, as in Figure [IJ which does not 
match the outer edge of the nonlinear tail in our fiducial 
model. Instead, the latter is well matched by rotating the 
former by 75° in the counterclockwise direction, shown 
as the dashed curve. 

4.2. Quasi- steady Density Wakes 

Figure 0] shows the quasi-steady distributions of the 
density as well as the velocity field, seen in the inertial 
frame, on the x/R p = 1, y = 0, and z = planes for 
our fiducial model at t/t or b = 2. The colorbar labels 
the perturbed density in logarithmic scale. The size of 
the short line segment shown in the upper left corner 



of each panel, corresponding to twice the sound speed, 
measures the amplitudes of the velocity vectors. The 
density wake consists of a spherical envelope surrounding 
the perturber, a detached bow shock, and an extended 
low-density region at the rear side. The envelope is al- 
most hydrostatic, as evidenced by low-amplitude velocity 
vectors in the y — plane to which the perturber mo- 
tion is perpendicular instantaneously. The bow shock 
standing outside the orbit, in the z — plane, is just a 
curved version of that in the linear trajectory counter- 
part and extends to large radii R = (x 2 + y 2 ) 1 ^ 2 in the 
orbital plane. On the other hand, the (almost planner) 
shock located inside the orbit gradually weakens at small 
R and terminates at R/R p ~ 0.5-0.6 since the speed of 
the density wake corotating with the perturber becomes 
subsonic there. The low-density regions delineating the 
tail seen at t/t OI b = 1 (Fig.[3)f) are less apparent in Fig- 
ure 2] as they diffuse out over time by pressure gradients. 
While the velocity just ahead of the perturber is parallel 
to the direction of the instantaneous perturber motion, 
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Fig. 6. — Density profiles of steady-state wakes as functions of the distance I from the perturber center along the azimuthal direction 
for (a) models with fixed TZ = 0.05 and varying B and (b) models with fixed B = 1 and varying 1Z. All the models have M = 2. In (a), 
the dotted lines give the respective density distributions under the assumption of hydrostatic equilibrium, while the numerical results are 
drawn as the solid lines. Note that varying TZ < 0.1 changes the density only in the very central parts of the hydrostatic envelopes. 



the direction of gas motions at the perturber location 
and the immediate behind it are inclined about 30°-50° 
from the tangential direction. This implies that the gas 
in the downstream side retains the information on the 
direction of perturber motion in the past. 

Figure [5] displays changes in the quasi-steady density 
wakes with varying M on the orbital plane for models 
with K = 0.02 and B = 1 (A = 20). In each panel, 
the perturber located at x/R p = 1 and y = is mov- 
ing along the white circle in the counterclockwise direc- 
tion. When M = 0.5, the density distribution is circular 
near the perturber and becomes eccentric away from it, 
with the minor axis parallel to the line of motion. Still, 
there exists a slight excess of the perturbed density in 
the rear side, giving rise to a non- vanishing (but small) 
DF force. Compared to the M. — 0.5 model, the model 
with A4 =0.8 has a larger density excess at the back- 
side. While the perturber is subsonic in this model, a 
strong gravitational pull is able to accelerate the back- 
ground gas to a near-transonic speed, creating a weak 
shock outside the orbit. The shape of the density wake 
in the M. = 1.2 model is similar to that for the M. = 2.0 
model, although the former has a larger peak density and 
a larger distance of the detached bow shock. As M. in- 
creases further, both the detached shock distance and the 
size of the hydrostatic envelope decrease. The ring-like 
low-density gap at R/R p ~ 2 in Figure EH is a remnant 
of the counterstream that was stretched along the orbit. 
With large values of M. and B, the counterstream in this 
model was so strong that the deep gap has not been hllcd 
in yet. 

In order to check whether the density wake near the 
perturber is in hydrostatic equilibrium, we plot in Figure 
[6] the density profiles of steady-state wakes as functions 
of distance I measured from the perturber location along 
the azimuthal direction in models with M. — 2. The 



solid lines in Figure [6]a give the simulation results for 
perturbers with fixed 1Z = 0.05 and differing B. The 
sharp discontinuity in each curve is due to the presence 
of the bow shock. The dotted lines draw the density 
distributions of the hydrostatic spheres 
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where p c and a c denote the density and adiabatic sound 
speed at I = (i.e., the perturber center), respectively 
(KK09). These are in good agreement with the den- 
sity profiles of the steady-state wakes in our simulations, 
suggesting that the envelopes are indeed almost hydro- 
static. Figure shows how the density profile varies 
with 1Z. When B = 1, the nonlinear effect is not signif- 
icant in models with 1Z > 0.33 since they have r] A < 1; 
the density jump occurring at I ~ in the 1Z = 0.4 
model is more like a Mach wave rather than a shock. As 
1Z decreases (or r\ A increases), the density wake becomes 
increasingly nonlinear, producing a hydrostatic envelope 
and a detached bow shock. Note that the shock loca- 
tion is unchanged as long as 1Z < 0.2 (or t) a > 2). When 
the density is highly nonlinear, the strong gravity as- 
sociated with small r s increases the density only in the 
small regions within < (1 — 2)r s . Since this core region 
is spherically symmetric, the resultant DF force becomes 
independent of 1Z provided r\ A > 2. We will show this 
directly in i )4. 41 below. 

Unlike in the linear-trajectory cases where a perturber 
travels through an undisturbed medium, it is inevitable 
that a circular-orbit perturber enters the backside of its 
own wake before completing one orbit. In the circular- 
orbit cases, therefore, the density ahead of the bow 
shock differs from the undisturbed density po and de- 
pends considerably on B, as Figure [6]a illustrates. The 
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Fig. 7. — Modified background density p\ ahead of the bow 
shock in models with 1Z = 0.05. Various symbols and errorbars 
give the means and standard deviations of the preshock density in 
the regions with Al/R p = 0.5—1.5 from the d eta ched shock. The 
solid line plots the best fit given by equation 1101 . 

enhanced preshock density can alternatively be inter- 
preted as arising from the response of the background 
gas to the orbit-averaged gravitational potential of the 
perturber (<J> P ) = J*° Ib <& P {t)dt/t or \>, which is centered 
at the orbit center. In this sense, the enhanced preshock 
density pi acts as a modified background density to a 
circular-orbit perturber with the gravitational potential 
$ P — ($ P ). To quantify p\, we take a spatial average of 
the preshock density in quasi-steady state over the re- 
gion with (Z — S)/Rp — 0.5-1.5 along the orbit, where 
6 denotes the location of the detached shock. Figure [7] 
plots the resulting p\ for models with 1Z = 0.05, with er- 
rorbars representing the standard deviations. The solid 
line is the best fit 
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showing that the change in the background density to 
circular-orbit perturbers is almost linearly proportional 
to B. Since pi is quite insensitive to M, the M- 
dependence of pi in equation (fTU|) can be ignored. 

4.3. Detached Shock Distance 

In laboratory experiments and hydrodynamic theo- 
ries, supersonic flows over a blunt-body generate a de- 
tached bow shock when the maximum angle allowed in 
the postshock flows is smaller than th e nose angle (e.g., 
iLiepmann fc Roshkol fl957t IShul fl99l) . Although a per- 
turber in our models does not hold any defined surface 
and merely provides gravitational perturbations to the 
background g spherical envelope formed around it 
nevertheless sends off pressure waves into the upstream 
direction, acting as an obstacle to the incoming gas. 
Since the nose-angle of a spherical body is 90° , the shock 
should be detached in our nonlinear supersonic models. 
Figure plots the time evolution of the detached shock 
distance 5 for some selected models with 1Z = 0.02 and 
M. = 2. Note that 6 initially increases rapidly and then 
saturates at t/t or h ~ 0.5 to an equilibrium value with 



some fluctuations. The quasi-steady value of 6 becomes 
larger as the perturber mass increases. 

For massive perturbers moving on a straight-line tra- 
jectory, KK09 found that the detached shock distance is 
given empirically by 



Sno = 



GM r , 



a 2 (M*-l) 



(11) 



This can be interpreted by the balance between the post- 
shock thermal energy (oc a 2 M 2 ) and the gravitational 
potential energy (oc —GM p /S) at the shock location for 
strong shocks (M ^ 1). To study how the circular obit 
affects the standoff distance of the shock, we take a time 
average of 6{t) over t/i rb ~ 1~3 and plot 6/6^ in Fig- 
ure [9] as various symbols. Again, errorbars indicate the 
standard deviations of 6(t) during this time interval. The 
numerical results show that the detached shock distances 
in the circular orbits are in general smaller than those in 
the linear-trajectory cases. This is likely because the con- 
tinuous change in the direction of perturber motion ef- 
fectively reduces the forward momentum of the incident 
gas near the perturber, as illustrated by the spatially- 
varying velocity field in the wake shown in Figure 2J In 
this case, the ram pressure of the incoming flow is able 
to push the shock front toward the perturber, making 
6 smaller. The modified background density and sound 
speed discussed in ^4. 21 might affect 6 as well. We fit our 
results using 



1 - 0.8A, 
0.14 + 0.2(1 



logA) : 



A < 0.1, 
0.1 < A < 10, 



(12) 



with A = 7 ]B /(M 2 - l) - 3 
lines. As expected, 6/6, 



which is plotted as dotted 
— > 1 as A — > corresponding 



to R p — > oo, although models with A < 0.1 are weakly 
nonlinear and have 6 comparable to or less than r s . Since 
6oo/Rp = '^B^ equation (|12p indicates that 5 increases 
with A, although its increasing rate slows down as A — > 
10 where 6 becomes comparable to R p . 

4.4. Drag Force 

For our nonlinear models, we calculate the DF force on 
a perturber due to its own induced wake. Figure [8^ il- 
lustrates the temporal changes of the dimcnsionlcss drag 
force X = F/[4irp {GM p /a ) 2 } in models with M = 2 
and 1Z = 0.02 and varying B. The drag force reaches 
a quasi-steady value in less than one orbit, although it 
exhibits some fluctuations for large B. This is markedly 
different from the cases of straight-line trajectories where 
the drag force increases logarithmically with time. The 
quasi-steady value of X is almost independent of B when 
B < 0.05, decreases with B for 0.05 < B < 2, and becomes 
again insensitive to B when B > 2. The behavior of X 
with B is a consequence of the circular orbit combined 
with the nonlinear effect. The size of a hydrostatic en- 
velope in a nonlinear density wake, as measured by the 
detached shock distance, becomes larger with increas- 
ing B, which makes the drag force smaller owing to the 
front-back symmetry in the wake. At the same time, the 
circular orbit draws the background material toward the 
orbit center and increases the density in the preshock 
region, which in turn tends to increase the drag force. 
For models with B < 2 (so that pi/po $ 2), the nonlinear 
effect dominates and X decreases with B. When B > 2, 
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Fig. 8. — Temporal evolution of (a) the detached shock distances <5 and (6) the DF forces F for models 1Z = 0.05 and M = 2. Both 5 and 
F increase rapidly with time initially and then saturate at t/t OI b ~ 0.5 to equilibrium values with some fluctuations. The dimensionless 
drag force decreases with B for 0.05 < B < 2 due to the nonlinear effect. 

instead of po. The dotted line marked with rj A = 2 de- 
marcates the parameter space into two parts such that 
the models shown in the lower-left region are fully non- 
linear with rj A > 2, while those in the upper- right region 
are in the linear or weakly nonlinear regime. The nor- 
malized drag force on weakly-nonlinear perturbers with 
0.7 < r) A < 2 is usually larger than that on low- mass per- 
turbers with rj A < 0.7 because the density wakes in the 
former are distributed closer to the perturber (KK09). 
In highly-nonlinear models with r\ A > 2 (or r\ e > 27?.), 
on the other hand, the DF force becomes essentially in- 
dependent of 1Z. This is of course because the change 
in r s affects the wake only within a small region close to 
the perturber where the density is spherically symmet- 
ric, making a negligible contribution to the drag force 
(see aOl . 

Figure [TUb plots F/F\i n versus r] A for the same models 
shown in Figure [TUb . Note again that po in equation (J]) 
for Fi; n is replaced by pi. For given B (so that increas- 
ing r\ A corresponds to decreasing 1Z), F/Fu n decreases 
logarithmically with increasing r\ A > 2. This is because 
r\ A oc IZ^ 1 for fixed B and M, leading to Fn n oc In r\ A , 
while F is constant for r\ A > 2. For fixed 1Z (so that in- 
creasing rj A corresponds to increasing B), F/Fy m oc rj A ' 5 , 
as illustrated by two dotted lines connecting the numer- 
ical results for 1Z = 0.05 and 0.025, respectively. Figure 
[TTI a plots F/Fy ln for all of our supersonic models with 
1Z = 0.05. Obviously, F/F\\ n « 1 for low-mass perturbers 
with rj A < 0.7. The dotted line is the best fit 

F 

F\in 

to the numerical results for nonlinear, circular-orbit per- 
turbers, which is very similar to equation for massive, 
linear-trajectory perturbers. 

In order to apply equation (|13p to estimate the de- 
cay timescale, one needs to specify the softening radius 




Fig. 9. — Quasi-steady values of the detached shock distances 5 of 
circular-orbit perturbers relative to those, 5^ , of linear-trajectory 
perturbers. Various symbols and error bars give the means and 
standard deviations of 8/ < 5qq over t/t or b = 1~3. The dotted curves, 
expressed as equation 1121 . are the best fits to the numerical results. 



on the other hand, the enhancement of the background 
density becomes significant enough to compensate for the 
nonlinear effect, making X unchanged with B. As we will 
show below, the DF force scales well with B when pi/po 
is taken into account. 

Compared to the linear-trajectory cases, the circular 
orbit introduces an additional length scale, R p , so that 
it is interesting to explore the dependence of the drag 
force upon 1Z. Figure 110b plots as various symbols and 
errorbars the mean values and standard deviations of the 
drag forces averaged over t/i rb — 2-4 for models with 
A4 = 2. Note that the drag force is normalized by p\ 
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Fig. 10. — Normalized DF forces averaged over t/t or b = 2—4 for 
models with M = 2 and various B and 7?. as functions of (a) 7?. and 
(b) t\a. The dotted line corresponding to t\a = 2 in (a) separates 
the domain into two parts such that the models at the lower-left 
part are fully nonlinear. Dotted lines in (b) have a slope of —0.5 
and connect the numerical results for 1Z = 0.05 and 0.025 models 
quite well. 

which is not directly observable in most astronomical sys- 
tems. In addition, equation (|13|) is applicable only when 
1Z = 0.05 since F on highly massive perturbers is inde- 
pendent of r s , while Fy ln diverges as r s — > 0. For prac- 
tical use in various situations, therefore, it is desirable 
to have an expression of F that does not rely on r s . For 
this purpose, we plot in Figure ITTb the dimensionless DF 
forces as a function of r\ B for all supersonic models with 
7Z = 0.05. The dotted line is our fit to the nonlinear part 
using 

4ir Pl {GM p ) 2 



F = 



V 2 

p 



7 

— ), for Vb > 0.1, 



(14) 



with pi given by equation (|10[) . For the models with 



r\ B < 0.1, equations ([T]) and (J5J with V p t 
a good estimate of F that depends on r s 



2R„ provide 



5. SUMMARY AND DISCUSSION 

DF due a gaseous medium may play a central role in 
removing angular momentum from astronomical objects 
in orbital motions, causing them to spiral in toward the 
orbit center. In this paper, we use three-dimensional 
hydrodynamic simulations to explore the gravitational 
wake and the associated drag force on a very massive 
perturber with mass M p moving at speed V p on a circu- 
lar orbit with radius R p . This work extends our pre- 
vious studies that considered low-mass, circular-orbit 
perturbers (KK07) and high-mass, linear-trajectory per- 
turbers (KK09). The background medium is assumed 



Fig. 11. — (a) DF forces on supersonic perturbers with 7Z = 0.05 
normalized by the linear drag force as a function of t\a. Various 
symbols and their errorbars give the means and standard deviations 
of F/Fn n for t/t or b = 2-4. The dotted line indicates a slope of 
—0.5. (b) The same data as in (a) but with different normalization 
and as a function of r\s ■ The dotted line shows the best fit to the 
numerical results for t)b > 0.1 and is given by equation 114> . 

to be uniform with density po and the sound speed an, 
adiabatic with index 7 = 5/3, non-self-gravitating, non- 
rotating, and unmagnetized. The perturber is repre- 
sented by a Plummer sphere with softening radius r s . 
Our models are characterized completely by three di- 
mensionless parameters: TZ = r s /R p , M = V p /ao, and 
B = GM p / (a^Rp). We run a total of 65 models that 
differ in it, Ai, and B, and obtain the dependence of the 
drag force on these parameters. 

For our supersonic models, the density wake of a mas- 
sive perturber reaches a quasi-steady state typically in ~ 
0.5i or b, about 10 times faster than the linear-trajectory 
cases. A quasi-steady density wake consists of a spheri- 
cal envelope surrounding the perturber, a detached bow 
shock in the upstream direction, and an extended low- 
density region in the rear side. The envelope is almost 
hydrostatic, providing a negligible contribution to the 
DF force. Compared to the linear-trajectory cases, an 
orbit-averaged gravitational potential of a circular-orbit 
perturber is able to gather the background gas toward 
the orbit center, effectively increasing the background 
density ahead of the bow shock. Equation (fTO)) describes 
the modified preshock density p±, showing that density 
enhancement in the preshock region is almost linearly 
proportional to B and depends weakly on M. 

Strong perturbations sent off from a massive perturber 
develop into a bow shock through which the upstream su- 
personic flow becomes subsonic. Since the density wake 
in a quasi-steady state forms a pattern that is corotating 
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with the perturber at a constant angular speed, the in- 
cident flow outside the orbit is supersonic to a perturber 
with M. > 1. This naturally produces a bow shock out- 
side the orbit that is simply a curved version of that in 
the linear-trajectory counterpart. Formation of a weak 
bow shock outside the orbit is possible even for a sub- 
sonic perturber with M. > 0.8 for which the incident flow 
is accelerated to near-transonic speeds. On the other 
hand, the motion of the steady-state wake is subsonic to 
the gas near the orbit center, even when M. > 1. This 
makes the bow shock gradually weaken at small radii, 
eventually terminating somewhere inside the orbit. The 
detached shock distance measured along the tangential 
direction is generally smaller than in the liner-trajectory 
counterpart since the continuous change in the direction 
of the perturber motion effectively reduces the forward 
momentum of the wake; equation (|12[) gives the alge- 
braic fits to the numerical results for the detached shock 
distance. 

Unlike in the linear-trajectory cases where the drag 
force increases logarithmically with time, the DF force in 
the circular orbits approaches a quasi-steady value in less 
than one orbit. The quasi-steady drag force is essentially 
independent of 1Z as long as rj A > 2 (or r\ e > 2TZ) since 
the change in r s modifies the density wake only within 
~ (1 — 2)r s from the perturber center, where the wake is 
almost spherically symmetric and thus has a negligible 
contribution to the drag force. For sufficiently massive 
perturbers, the presence of a near-hydrostatic envelope 
in the density wake also makes the nonlinear drag force 
smaller than the linear estimate. Since the choice of r s 
is uncertain in many practical applications, we fit the 
nonlinear drag force on a circular-orbit perturber using 
equation (fT4")l that does not involve r s . 

That the nonlinear drag force in circular orbits is in- 
sensitive to the softening radius seems inconsistent with 
the results of KK09 for linear-trajectory perturbers. But, 
it does not. For the latter case, r s is in fact a lone length 
scale relative to which other length scales such as the 
Bondi radius and the distance traveled by the perturber 
are measured. It is then difficult to separate the effect of 
r s on the drag force from that of M p because the models 
depend on them only through the dimensionless param- 
eter A. Had we run simulations by using dimensional 
quantities in KK09, we would probably have gotten the 
results that F oc Mp 5 independent of r s for models with 
r) A ^> 1, similarly to the results of the present paper. 
Given that Fn a oc Mp ln(V p t/r s ), the time-averaged val- 
ues of F/Fun, when expressed in terms of r\ A in KK09, 
somehow picked up the power-law dependence on M p 
more clearly than the weak logarithmic dependence on 
r s . The circular orbit breaks the degeneracy between r s 
and M p by introducing another length scale, R p , allow- 
ing to separate the effect on the DF force of r s from that 
of M p . The similarity between equation (p~3|) obtained 
for fixed 1Z in the current work and the result of KK09 
suggests that the latter should be interpreted as the de- 
pendence of F on M p for fixed r s . 

Now we compare the results of our idealized mod- 
els with those from realistic simulations fo r the orbital 
decay of SMBHs at galaxy centers (e.g.. lEscala et al.l 
12004 120051: [Dotti et alil2006l 120071: [Maver et al.ll2007t see 
also IColpi fc Dottil 120091 and references therein) . While 



these authors showed that SMBHs inspiral rapidly in 
~ 10 6 — 10 7 yrs due to the gaseous DF force, they 
also found that the DF force in supersonic models is 
~ 1.5 times smaller, and depends on the black hole 
mass les s sensi tively, than the analytical predictions o f 
lOstrikerl (|1999h formula (e.g., lEscala et alJl2004 l2005| ). 
The discrepancies between the numerical and analytical 
results are most likely due to the nonlinea r effect. For 
instan ce, one of the models considered by lEscala et al.1 
(|2005f) has a black hole with mass M p = 5 x 10 7 M and 
softening radius r s = 4 pc moving initially with speed 
V p ~ 100 km s -1 on a circular orbit with R p = 200 pc 
through a medium with sound speed do ~ 70kms~ 1 , 
corresponding to B = 0.2, M = 1.4, and TZ = 0.02. 
For this model, equation (|14p compared with equation 
(HI) would give Fiin/ F ~ 2, consistent with the results of 
lEscala et al.l j2005). Since the orbital decay timescale 
is given by Tdecay = M p V p /F, equation (fT4"| predicts 
Tdccay oc M~ - 5 for sufficiently massive perturber, which 
is more consistent with Tdecay oc M~ 3 inferred from 
the results of lEscala et all (|2004 l2005f) than the linear 
prediction Td cca y oc M~ . The difference between the re- 
sults of the current paper and lEscala et alj (|2004 l2005h 
is presumably due to the effects of density stratification, 
rotation, self-gravity in the background medium, which 
are not considered in the present work. 

Another important issue regards the proper choice of 
the softening length of a point-mass perturber employed 
in numerical simulations for the gaseous DF. As men- 
tioned in fJTJ the limited numerical resolution commonly 
requires to take r s a few orders of magnitude larger than 
the realistic size. Are such large values of r s acceptable 
in calculating the orbital decay time accurately? Based 
on our numerical results, the answer is yes, provided 
r\ A > 2 (or rj e > 21Z). Otherwise, a large value of r s 
would cause the gravity near the black hole to be re- 
duced significantly. The resulting density wake would 
then remain in the linear regime, making the drag force 
overestimated considera bly. For the decay of SM BHs, 
the models considered in lEscala et al.l (|2004 l2005f) have 
rj A ~ 10, well above the required lower limit. Some mod- 
els with M v ~ 10^ - 10 7 M in iDotti et all (|2007l ) and 
iMaver et~a l. ( 2007) have r] A ~ 1, so that their choices of 
r s ~ 0.1 — 0.2 pc are marginally acceptable. For the evo- 
lution of a main-sequence comp anion with M = . 4 
in a common-envelope binary, Sandqui st et al.l (|1998l ) 



took 



2.3 x 10 cm, about 10 times larger than 



the real size. Since this corresponds to r) A ~ 20, our re- 
sults suggest that the softening radius did not affect the 
decay time in their numerical simulations. 

Finally, we remark on the several assumptions made 
in the current work. First, we consider a single per- 
turber moving on a circular orbit. If the resulting drag 
force (eq. [14]) is to be applied to the decay of double 
black holes, one has to take into account the drag force 
from the wake of its co mpanion located at the other side 
of the orbit, as well. iKim et aTl (|2008f) found that for 
equal-mass perturbers, the ratio of the DF force from 
the companion wake to that from its own induced wake is 
about ~ 10 — 50% for supersonic cases, depending on the 
Mach number, when the perturbers have so low masses 
that the wakes are in the linear regime. The nonlinear 
effect on massive perturbers has yet to be explored. Sec- 



14 



W.-T. Kim 



ond, the current study assumes an adiabatic equation of 
state with index 7 = 5/3. It is not well known what 
the most appropriate value of 7 is in the background 
medium, but it would be 7 = 5/3 if the gas is fully 
ionized or 7 — 7/5 for preferentially molecular gas, un- 
der the assumption that radiative heating and cooling 
are un important. The numerical results of iMaver et al.l 
( 2007) indicate that the orbital decay is more effective for 
a softer equation of state, suggesting that the DF force 
may depend sensitively on 7. Third, while we consider 
a static background medium with uniform density, it is 
more likely to be strat ified in real situations. Fo r a col- 
lisionless background, Uust fc Penarrubial ()2005l ) found 
that a density gradient induces an additional drag force 
in the lateral direction of the perturber motion, amount- 
ing to about 10% of the drag in the backward direction. 



We also have ignored the rotation, self-gravity, buoyancy, 
and turbulent motions of the background medium in the 
present paper. It is interesting to see what effect each 
of these physical ingredients makes, which will direct our 
future research. 
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